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Effect of electron diamagnetic drifts on cylindrical double-tearing modes 

Stephen Abbott^ and Kai Germaschewski^ 

University of New Hampshire, Durham, NH 03824 

Double-tearing modes (DTMs) have been proposed as a driver of ‘off-axis sawtooth’ crashes in reverse magnetic 
shear tokamak configurations. Recently differential rotation provided by equilibrium sheared flows has been 
shown capable of decoupling the two DTM resonant layers, slowing the growth the instability. In this work 
we instead supply this differential rotation using an electron diamagnetic drift, which emerges in the presence 
of an equilibrium pressure gradient and finite Larmor radius physics. Diamagnetic drifts have the additional 
benefit of stabilizing reconnection local to the two tearing layers. Conducting linear and nonlinear simulations 
with the extended MHD code MRC-3d, we consider an m = 2, n = 1 cylindrical double-tearing mode. We 
show that asymmetries between the resonant layers and the emergence of an ideal MHD instability cause the 
DTM evolution to be highly dependent on the location of the pressure gradient. By locating a strong drift 
near the outer, dominant resonant surface are we able to saturate the mode and preserve the annular current 
ring, suggesting that the appearance of DTM activity in advanced tokamaks depends strongly on the details 
of the plasma pressure profile. 


I. INTRODUCTION 

The promising properties of the reversed magnetic 
shear (RMS) tokamak configuration have led to a recent 
increased interest in the double tearing mode. RMS con¬ 
figurations result in safety factor profiles that are non¬ 
monotonic. If two rational surfaces of the same q exist 
near each other within the plasma, they may couple to¬ 
gether via ideal MHD scale processes to form a single 
double-tearing mode (DTM). Linearly the interaction of 
the two surfaces creates a self-driven reconnecting insta¬ 
bility that depends weakly on resistivity.^^ Nonlinearly the 
DTM can potentially disrupt the annular current ring of 
RMS devices,!^ generate strong sheared flows,and re¬ 
lease large bursts of kinetic energy.® As such, they are a 
proposed driver of off-axis sawtooth behavior. 

One means of stabilizing double-tearing mode activ¬ 
ity is the application of differential rotation to the two 
DTM layers. In slab Cartesian simulations equilibrium 
sheared flow has been shown to interfere with the cou¬ 
pling between the resonant layers to result in two de¬ 
coupled, drifting single tearing modes.® Further increase 
in the flow amplitude generates Alfven resonance layers 
that couple to the tearing surfaces, increasing or decreas¬ 
ing the mode growth depending on their proximity. Non- 
linearly these Alfven resonances may shield the plasma 
core and suppress DTM mode growth.®^ The appearance 
of such layers requires, however, flow amplitudes near 
the in-plane Alfven speed and shears near the threshold 
for Kelvin-Helmholtz instability,® potentially triggering 
greater instability. Thus we are motivated to explore al¬ 
ternate means of providing differential plasma rotation. 

Diamagnetic drifts emerge as a result of finite Lar¬ 
mor radius physics in the presence of a pressure gradient 
such as the internal transport barriers (ITBs) observed in 
RMS plasmas.® They have long been studied as means of 
stabilizing reconnecting modes, and have several advan¬ 
tages over equilibrium flow. In particular, diamagnetic 
drifts local to the reconnecting layer interfere with the 
conversion of magnetic energy to kinetic.® This local ef¬ 
fect has been shown to saturate the m = 1 kink-tearing 


mode in conventional tokamaks,!^ leaving finite sized is¬ 
lands during incomplete sawtooth crashes.^ 

The influence of both pressure gradients and diamag¬ 
netic drifts on double-tearing modes has been considered 
previously by other authors. In resistive, reduced MHD 
simulations Zhao et al.^ examined the impact of equi¬ 
librium pressure gradients on a cylindrical DTM with a 
small inter-resonant spacing. Their results show that the 
pressure gradient modifies the dependence of the DTM on 
resistivity, causing variations in the spectrum of modes at 
a given surface. In this work we will expand their study 
to more widely spaced modes and a wider variety of pres¬ 
sure gradients, as well as introduce finite Larmor radius 
effects. Maget et al.l^^ considered neoclassical effects on 
a DTM in toroidal simulations, targeting Tore Supra ex¬ 
periments. They found some mode numbers were sta¬ 
bilized by the addition of diamagnetic drifts, whereas 
others were enhanced due to toroidal effects. Their sim¬ 
ulations target, however, specific discharges and do not 
consider variations in the pressure and drift profiles, thus 
do not illuminate the role of differential rotation in DTM 
evolution. 

In this work, we use Hall MHD simulations to examine 
the impact of diamagnetic drifts on a cylindrical m = 2, 
n = 1 double-tearing mode, considering both the ability 
of an electron diamagnetic drift to decouple the two tear¬ 
ing layers and to stabilize them once separated. To this 
end, we structure this paper as follows. In Section |TT| we 
introduce the Hall MHD model and describe our simu¬ 
lation code MRC-3d. Section [m] defines the equilibrium 
safety factor and density profiles used for this study. In 
Section ||V] we report the results of linear resistive and 
Hall MHD simulations. We find that the addition of a 
pressure gradient to our equilibrium destabilizes an ideal 
MHD instability that competes with the stabilizing ef¬ 
fects of the diamagnetic drift. As a consequence, we are 
able to decrease the linear DTM growth rate only by lo¬ 
cating a strong diamagnetic drift at the dominant, outer 
rational surface. We use this result to choose charac¬ 
teristic profiles for nonlinear Hall MHD simulations in 
Section |Vj and show that the DTM may be saturated 
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before disruption of the annular current ring. Finally, 
we summarize our results in Section Ell and discuss the 
consequence of this work for advanced tokamaks. 


II. HALL MAGNETOHYDRODYNAMIC MODEL 

Our simulation code MRC-3d uses a standard Hall MHD 
model. 
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where p = (1 + T)pris the pressure and t = TijT^ is the 
ratio of the ion to electron temperatures. In this work 
we focus on the cold-ion regime (r = 0), thus excluding 
ion diamagnetic drifts. 

The simulation code MRC-sdSl implements the above 
model in a fully conservative, finite-volume scheme sim¬ 
ilar to ChacorpSl^ with additional di scale Hall and elec¬ 
tron pressure gradient terms. Lengths are normalized 
to the cylinder radius, magnetic fields to the asymp¬ 
totic in-plane magnitude, and velocities to the in-plane 
Alfven speed. All other normalizations follow from these. 
Data management and implicit time integration are ac¬ 
complished via the interface in the LIBMRC 

computational library.l^ 

We conduct simulations in 2D helically symmetric 
cylindrical geometry. Derivatives in radial 0 < r < 1 
and poloidal 0 < 0 < 27r coordinates are discretized di¬ 
rectly. The cylinder is assumed to be periodic in z with 
a length of 2ttR where i? = 10 is the major radius of an 
approximately equivalent torus with inverse aspect ratio 
e = 0.1. Derivatives in the axial coordinate (j) = zjR are 
taken to be dcj) = c~^d6 where t = n/m is the twist of the 
helix. Thus for analysis we define the helical coordinates 
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where f and u represent the two dimensional perpendic¬ 
ular plane and h is directed along the helix. The helical 
flux function tjj* is then defined as 

( 10 ) 


To reduce the computational costs of this work we use 
resistivities on the order of p ^ 10“®. This unrealistically 
large level of diffusion causes the equilibrium to decay on 
a time scale comparable to the growth time of resistive 
instabilities. MRC-3d features a mechanism to prevent 
equilibrium decay that is equivalent to the addition of a 
source electric field in Ohm’s law (Eqn. . In nonlinear 
simulations we enable this mechanism and prevent the re¬ 
sistive decay of the equilibrium. We have confirmed that 
the major results of this work persist when the source 
electric field is disabled, and will note explicitly when it 
impacts the mode behavior. 


III. EQUILIBRIUM 

To generate an equilibrium with two nearby q = 2 
rational surfaces we use the safety factor profile given by 
Bierwag^^ 

- N 1 /w(r) 

ro = rA\[m/inqo)r^^-^ - l|-i/[2»hH] 

w{r) = Wq + Wir'^ 

Fi{r) = 1 + fi exp {-[(r - rii)/ri 2 ]^} 
with the following constants: 

ta = 0.655, Wo = 3.8824, wi = 0 


/i = -0.238, rii = 0.4286, ri2 = 0.304 


( 12 ) 


We set qo = 2.5, resulting in two q = 2 rational surfaces 
located at 0 < rsi < rs 2 < 1, spaced a distance D = 
rs 2 — « 0.26 apart. 

Diamagnetic drifts require the introduction of a pres¬ 
sure gradient. We use a monotonic density profile to rep¬ 
resent an internal transport barrier as given by ZhacP^ 


p(r) = iVo 1 - (1 - iVb) 


tanh(ro/dAr) -|- tanh[(r — ro)/(5Ar] 
tanh(ro/djv) + tanh[(l — ro)(yAr] 

(13) 

We fix the core density with Nq = 1 and vary tq I 
and Nb to change the center, width, and edge magnitude 
of the density profile. For simplicity we take the equi¬ 
librium electron temperature to be a constant Tq = 1.0, 
and assume cold-ions r = T^/Te = 0 . We initialize 
the equilibrium magnetic field Bq by specifying a den¬ 
sity profile and iteratively refining Bgo and B^o toward 
the above safety factor profile, subject to the constraints 
of force-balance and that B^o ~ 10. This equilibrium 
has a plasma parameter of f3 « 0.01 so that /3 ~ e^, 
consistent with the standard tokamak ordering assump¬ 
tion. The equilibrium safety factor and example density 
gradient are shown in Figure 

This pressure gradient produces (in the cold ion regime 
with a constant electron temperature Tq = 1.0) an elec¬ 
tron diamagnetic drift given by: 

mdiBh dp 


a;*(r) = k ■ v^e = —■ 


rpB"^ dr 


B = •Kh + Bhh 


(14) 
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FIG. 1. Equilibrium safety factor profile (solid, left axis) with 
two nearby q = 2 rational surfaces, indicated by dot-dashed 
vertical lines, and an example pressure profile (dashed, right 
axis). 

where we have assumed the perturbation wave vector is 
k = m/rO — n/Rz, commensurate with our helical sym¬ 
metry. 

We choose three classes of density prohles, shown in 
Figure which center the maximum gradient at differ¬ 
ent locations ■Csi < rg < rs 2 - Equal drift profiles have 
the peak gradient centered between the two q = 2 res¬ 
onant surfaces at rg = (rsi -|- rs 2 )l 2 and a broad width 
of 6n = 0 . 2 , producing equal w* at both singular lay¬ 
ers. Inner drift profiles have a narrow pressure profile of 
6m = 0.05 centered at rg = r^i, providing a strong drift at 
the inner rational surface and negligible w* at the outer. 
Finally, outer drift prohles are localized near rg = rs 2 
with 5m = 0.05 so that the inner rational surface expe¬ 
riences negligible w*. Equal drift prohles demonstrate 
the stabilizing effects of local diamagnetic drift on both 
surfaces. The inner and outer drift prohles produce a 
differential diamagnetic drift Aw* = |w*(rsi) — w»(rs 2 )|, 
resulting in an additional differential rotation effect and 
asymmetric local stabilization of the two surfaces. The 
impact of diamagnetic drifts in more realistic ITB-like 
prohles will likely be some intermediate form of these 
three prototypical equilibrium types. 

IV. LINEAR DRIFTS 

The basic consequences of introducing diamagnetic 
drifts to DTM evolution are clearly observable during 
the linear phase. MRC-3d includes a one dimensional, lin¬ 
earized form of the model given in Section |llj which we 
will use for this portion of the study. We treat deriva¬ 
tives in r using hnite volume discretization and apply 
the Fourier ansatz F{r,d,z) = f (r) exp {m9/r — nz/R) 
to derivatives in 9 and z, where the poloidal and toroidal 
mode numbers m = 2 and n = 1 are chosen to capture 
the lowest (and fastest growing) harmonic. From initial 


value simulations we ht the growth rate "fn of the mode 
from the time evolution of the magnetic and momentum 
held amplitudes. To extract the mode drift frequency jj 
we apply a discrete Fourier transformation (DFT) to the 
time series output of the helical hux function fj*. 

The reduction to a one dimensional linear model allows 
us to easily conduct scaling studies of DTM behavior in 
the three equilibria types given above. For each prohle 
type, we run simulations over a range of diamagnetic drift 
values 0 < w* < 0.16. For each simulations we specify the 
center of the gradient rg and the desired drifts at the in¬ 
ner (a;*(rsi)) and outer (a;*(rs 2 )) q = 2 rational surfaces, 
then iteratively rehne the density height W and the mag¬ 
netic helds to produce the desired prohle. We then seed 
a small m = 2, n = 1 perturbation onto this equilibrium 
and run the initial value simulation using a resistivity 
r] = 1 X 10“®. To enhance numerical stability we add in 
small amounts of viscosity, particle dihusivity, and tem¬ 
perature diffusivity v = D = DT = 10“^ x 77 . These 
extra dissipation coefficients smooth noise in the linear 
simulations and allow easier analysis; we have conhrmed 
that they have a negligible impact on the measured linear 
growth rates. 


A. Ideal MHD instability 

Before examining the impact of the diamagnetic drift 
we hrst set di = 0 and consider the addition of a pressure 
gradient in resistive MHD. In Figure]^ we have plotted 
the dependence of the linear growth rate on the max¬ 
imum of the pressure gradient drP\r=ro- Although the 
diamagnetic drift is not present in these resistive simu¬ 
lations, we will continue referring to the three types of 
profiles as equal (rg = (rgi + rs2)/2), inner (rg = r^i), 
and outer (rg = rs 2 ) ‘drift’ configurations. 

Increasing the pressure gradient increases the growth 
rate for all three classes of profile, although most dra¬ 
matically for the ‘equal drift’ case. This dependence of 
the growth rate on pressure gradient suggests the m = 2 , 
n = 1 DTM couples to an ideal MHD instability, similar 
to the interaction between the m = I, n = 1 kink and 
tearing modes.^^ We verify the presence of this ideal in¬ 
stability by running a scaling study of growth rate with 
resistivity in three sample equilibria. Figure shows that 
in the presence of a pressure gradient there is a mini¬ 
mum growth rate below which 7 no longer varies with 
resistivity. This minimum value increases with increas¬ 
ing pressure gradient, and is not observed in the force-free 
equilibrium. While we have not conmleted the analysis 
of this ideal mode, Pritchett et al.^ showed in Carte¬ 
sian geometry that the DTM tearing layers couple to a 
slab-kink mode, the stability of which determines the de¬ 
pendence of the growth rate on resistivity. We propose 
that in cylindrical geometry the addition of a pressure 
gradient may cause this kink mode to become unstable, 
thus further driving the DTM growth. 

We note that Zhao et found a similar increase in 
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FIG. 2. Examples of the three pressure prohle types (solid lines) and the electron diamagnetic drifts they produce (dashed 
lines) from left to right: equal drift; inner drift; and outer drift. Vertical dash-dot lines indicate the locations of the q = 2 
rational surfaces. All three profiles results in a drift of cu* = 0.8 at one or both of the rational surfaces, respectively. 
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FIG. 3. Variation in linear m = 2, n = 1 DTM growth 
rates yrj with pressure gradient in resistive MHD. For better 
comparison with Figurej^we represent the pressure gradients 
using tUres, which is the electron diamagnetic drift each profile 
would produce at the specihed resonant surfaces were di =0.1 
rather than 0. The three classes of profile described in Sec¬ 
tion III are represented as: triangles-ro = {rsi+rs2)/2, (equal 
drift); squares-ro = r^i (inner drift); circles-ro = rs 2 (outer 
drift). 


the DTM growth rate in the presence of a nontrivial pres¬ 
sure profile. In their simulations, however, the pressure 
gradient increased the dependence of the m = 2, n = 1 
DTM on resistivity (7 oc The inter-resonance dis¬ 

tance oi D = 0.26 examined here is much larger than 
than the D = 0.06 mode considered by Zhao, suggesting 
that the impact of the pressure gradient may depend on 
the spacing between the rational surfaces. 


B. Diamagnetic drift effects 

Having established via resistive simulations that the 
pressure gradient introduces an ideal MHD instability, 
we now introduce finite Larmor radius effects and con¬ 
sider how the diamagnetic drift impacts the cylindri¬ 
cal DTM. We fix the ion inertial length at di = 0.1 in 
Eqn. which results in an ion-sound Larmor radius of 


-A - Force-free (3,f)(r=rQ) =0) 

(},p('-='-q)=-1-12 

(J,.p(r=ro)=-1.79 



FIG. 4. Dependence of the m = 2, n = 1 DTM on resis¬ 
tivity rj in the presence of a pressure gradient centered at 
fo = (f'si + rs 2)/2 (equal drift profile) for different values of 
the peak pressure gradient. Triangles drP|r=ro = 0; squares 
drPlr=ro = —1.12; CUcleS drPIr^ro = “1.79. 


Ps = y/Pdi « 0.014. This large ion scale is required 
to provide sufficient scale separation given our use of 
a large resistivity to enhance numerical stability. With 
di fixed, increasing the maximum pressure gradient in¬ 
creases the diamagnetic drift frequency. In Figure we 
show the effect of increasing diamagnetic drift at the in¬ 
ner (tq = Tsi), outer (rg = rs 2 ), or both q = 2 rational 
surfaces. 

Both the equal drift and inner drift profiles are dom¬ 
inated by the ideal MHD behavior observed in resistive 
MHD simulations. For w* < 0.05, equal drift at both 
q = 2 rational surfaces counterbalances the enhanced 
driving energy and the growth rate remains almost con¬ 
stant with increasing pressure gradient. The diamagnetic 
drift is, however, unable overcome the ideal mode at large 
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FIG. 5. Variation in growth rate 7 ij and mode drift frequency 
7 / with a diamagnetic drift of cn* at: both q = 2 rational sur¬ 
faces (triangles, equal drift); the rai surface (squares, inner 
drift); the rs 2 surface (circles, outer drift). Simulations are 
conducted in Hall MHD, di — 0.1, ps = 0.014, with pressure 
gradient as described in Sec. |III| 


pressure gradients and 7 ^ again tracks the resistive simu¬ 
lations. We do not observe any region of constant growth 
rate for the inner drift equilibria, and the DTM behavior 
is dominated by the ideal MHD driving for all pressure 
gradients. 

Localized w* at the outer resonant surface has a much 
stronger stabilizing affect on the DTM. At small drifts 
(w* < 0.025) the growth rate decreases slowly with in¬ 
creasing An inflection point is evident in the scaling 
near w* = 0.025, after which the growth rate decreases 
more rapidly and nearly linearly. The eigenmode at a 
drift of w* = 0.02 (ro = rs 2 , Nb = 0.949, Sn = 0.05) 
shows a shearing of the perturbed helical flux tjj* between 
the inner and outer q = 2 surfaces that is absent in re¬ 
sistive simulations of the same equilibrium (Figure |^. 
Similar shearing of the eigenmode has been observed as 
consequence of equilibrium sheared flow.l^ 

Differential diamagnetic drifts above the critical value 
of Aw* = w*(rs 2 ) = 0.025 cause decoupling of the recon¬ 
necting layers, i.e. the system acts predominantly as two 
independent, drifting, single tearing modes rather than 
a single double-tearing mode. Comparing two different 
times of an w* = 0.1 outer-drift simulation in Figure 
shows independent movement of the structure around the 
inner and outer rational surfaces. This decoupling is, 
again, similar to that observed in resistive MHD sheared 
flow studies.l^ The continued decrease of the growth rate 
above w* « 0.25 is not present in sheared flow equilibria; 
it is instead due to stabilizing effect s of the diamagnetic 
drift local to the singular layerThus the outer drift 
equilibrium manifests both the decoupling properties of 
equilibrium sheared flow and the reconnection inhibiting 
benefits of the diamagnetic drift. 

The scaling behavior in Figure shows significant dif¬ 
ferences depending on where we apply the diamagnetic 
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FIG. 6. Linear eigenmodes of the helical flux function tp* in 
the presence of pressure gradient centered at ro = rs 2 with 
Nb = 0.949. In restive MHD {di — 0, left) there is no diamag¬ 
netic drift; in Hall MHD {di = 0.1, right) the outer surface 
experiences a drift of w* = 0.02 while the inner surface does 
not, resulting in a shearing of the eigenmode. 
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FIG. 7. Two simulation times for a linear m = 2, n = 1 
DTM with pressure profile centered at ro = rs 2 producing 
a localized drift of w* = 0.1. The perturbed helical flux ip* 
at the later time (right) is not a simple rotation of the ear¬ 
lier time (left). Instead the inner and outer rational surfaces 
evolve independently, overlapping between the surfaces. The 
two tearing layers are decoupled from each other by the dif¬ 
ferential diamagnetic drift. 


drift. These differences are due to the inherent asymme¬ 
try between the two rational surfaces of the cylindrical 
DTM. Consider the growth rate 'ju and eigenmode drift 
frequency 7 / for the stabilized outer drift (ro = rs 2 ) and 
weakly destabilized inner drift (rg = rsi) profiles. Both 
feature a strong diamagnetic drift localized at one of the 
rational surfaces, producing a differential Aw* that re¬ 
sults in decoupling. Increasing w* at rg = Tsi does not, 
however, produce any measurable eigenmode drift ( 7 /). 
Even in the outer drift (rp = rs 2 ) case, where the pertur¬ 
bation does rotate, only one drifting eigenmode can be 
found via Fourier analysis. This behavior is in contrast to 
slab-Cartesian sheared flow studies where two oppositely 
drifting eigenmodes are observed post-decoupling.l^ 

The magnetic shear is much greater at the outer ratio¬ 
nal surface than the inner (see Fig. [^, and as a conse¬ 
quence the driving energy local to the rs 2 tearing layer is 
much larger. When the two layers are coupled they grow 
as a single mode, but the equilibrium asymmetry causes 
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the eigenfunction to be biased toward the outer rational 
surface (Figure [^. This surface largely dominates the 
DTM growth. When decoupled the single layer at rs 2 is 
the fastest growing mode; the slower mode at the inner 
rational surface cannot easily be detected in our initial 
value simulations. 

A diamagnetic drift localized near ro = r^i stabilizes 
and decouples the weaker inner rational surface. The 
outer, dominant surface does not experience any drift, 
and thus the inner drift pressure profile does not result 
in a measurable eigenmode drift ( 7 /). Nor does the domi¬ 
nant, outer surface encounter any of the stabilizing effects 
of the diamagnetic drift. As a consequence, ■jr is largely 
controlled by the destabilization of the ideal MHD mode. 

When the drift is instead localized near rg = rs 2 , the 
outer surface does rotate and we measure a finite 7 /. 
The dominant surface now experiences the stabilizing 
diamagnetic effects and the growth rate decreases. For 
w* > 0.1 the growth rate again slowly increases with in¬ 
creasing pressure gradient, and the measured mode drift 
frequency 7 / suddenly drops to a much lower value. Con¬ 
sidering Figure]^ the eigenfunction near the inner (rsi) 
singular layer is clearly visible. The local diamagnetic 
drift at rs 2 has thus stabilized the outer layer sufficiently 
that it now has a slower growth rate than the inner, un¬ 
stabilized tearing mode. Increasing the pressure gradi¬ 
ent beyond this value will have not further decrease the 
growth rate, as now the ideal MHD driving of the inner 
surface has become dominant. 

Our linear simulation results highlight two qualities of 
this m = 2 , n = 1 double-tearing mode which limit the 
stabilizing properties of the diamagnetic drift. Firstly, 
the DTM is strongly driven by the interaction between 
the q = 2 rational surfaces. Unless a profile provides 
some differential rotation effect to decouple the two tear¬ 
ing layers, the diamagnetic drift is not sufficient to over¬ 
come the ideal MHD driving of the increased pressure 
gradient. Secondly, the asymmetric magnetic shear in¬ 
herent to cylindrical geometry causes one of the q = 2 
rational surfaces to be dominant. For the eigenmode 
growth rate to decrease, this fastest growing surface must 
experience the stabilizing diamagnetic drift. 


V. NONLINEAR DIAMAGNETIC DRIFTS 

Based on the linear results of the previous section, we 
choose the drift of w* = 0.1 for nonlinear simulations. 
When localized near the outer rational surface this dia¬ 
magnetic drift resulted in the lowest observed growth 
rate. To better understand the nonlinear evolution of 
the DTM we will compare the w* = 0.1 outer drift pro¬ 
file to both the force-free DTM and equal and inner drift 
profiles with the same w*. 

The nonlinear evolution of double-tearing modes is 
commonly classified by the growth of the kinetic and per¬ 


turbed magnetic energies: 

Ek = J^pV-UdV (15) 

E^ = - J ^{SB- 6B + 25B ■ Bg) dV (16) 

The introduction of diamagnetic drifts in this work will 
cause the perturbations to rotate, thus the kinetic en¬ 
ergy of all the drifting systems will typically be larger 
than the force-free case. As a consequence, the absolute 
magnitude of is not as good a representation of DTM 
stability as Em- The general features of the kinetic en¬ 
ergy will, however, provide an indicator of the stage of 
DTM evolution. 

The kinetic and magnetic energy growth of the force- 
free, m = 2, n = 1 baseline is shown in Figure]^ The long 
period of nearly exponential energy growth represents the 
development of finite sized magnetic islands such as those 
in Figure In our simulations these magnetic islands 
do not develop the magnetic structure necessary for the 
explosive growth phase observed in higher mode num¬ 
ber DTMsp^ the kinetic energy of the force-free DTM 
in Figure approaches a maximum value at simulation 
time t « 265 smoothly. This maximum E^ corresponds 
to the separatrix at the inner rational surface merging 
with that of the outer, as show in Figure At this time 
the flux between the rational surfaces has been consumed 
by the reconnecting layers and the annular current ring 
is disrupted. Continued evolution beyond this point re¬ 
sults in the magnetic islands reconnecting completely and 
relaxation of the system. 

A particular feature of this moderately spaced, low 
mode number DTM is that the flux surrounding the mag¬ 
netic axis is consumed by reconnection at the r^i rational 
surface, causing the inner current sheets to merge across 
r = 0. Other simulations of off-axis sawtooth activity 
in TFTIlPhave shown similar behavior, however it is not 
typically obser ved in higher mode number or more closely 
spaced DTMs.l^^^ In this work we treat the separatrix 
merging event as a complete loss of system stability, and 
thus will not consider whether such highly symmetric be¬ 
havior is relevant to realistic devices. 

To evaluate the impact of electron diamagnetic drifts 
nonlinearly we consider examples from each of the three 
classes of pressure profile (equal, inner, and outer drift) 
that have an equilibrium drift of w* = 0.1 at both ra¬ 
tional surfaces, just the inner, or just the outer. To 
decrease simulation time we seed our 2D helically sym¬ 
metric simulations with an m = 2 , n = 1 perturbation 
approximately 1 x 10 “^ times the equilibrium field, and 
evolve the system using an implicit, second order Crank- 
Nicolson method. For these simulations we use a grid of 
2048 cells in r and 512 in 9, with a nonuniform distribu¬ 
tion in the radial coordinate that increases resolution at 
and within the rs 2 surface while decreasing it towards the 
r = 1 conducting wall boundary. To aid numerical sta¬ 
bility we choose a resistivity of ry = 2 x 10 “® and set all 
other dissipation coefficients to 1 x 10“®. We have sue- 
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FIG. 8. Nonlinear kinetic and energy growth of the DTM in 
the presence of various pressure profiles. The force-free profile 
has no equilibrium diamagnetic drift. The remaining profiles 
have drifts of oj* = 0.1: rg = (r^i + rgg)!^ - equal drift at 
both q — 2 surfaces; ro = r^i - localized at the inner surface; 
ro = rs 2 - localized at the outer surface. 


t = 180 t=26o 




FIG. 9. Contours of the helical flux function ip* show the is¬ 
land growth regime (left) and separatrix merging event (right) 
of the nonlinear, force-free DTM. 


cessfully convergence tested the following results in both 
spatial and time resolution, and also verified that the ex¬ 
tra dissipation coefficients do not significantly impact the 
mode behavior. 


A. Inner drift 


During the linear phase, locating a pressure gradient at 
the inner rational surface produced a marginal increase 
in the growth rate but otherwise did not significantly 
impact the mode evolution. Nonlinearly we find that 
the dominance of the rs 2 rational surface persists and a 
strong pressure gradient (and thus drift) at results in 
slightly faster growth of the magnetic perturbation when 
compared to the force-free system (Figure [^. We find 
that the island at the outer rational surface quickly grows 
large enough to interlock the layer at the inner surface, 
which has been previously been observed in nonlinear 



In-plane flow 



-0.50 0.00 0.50 


FIG. 10. Applying a diamagnetic drift of tu* =0.1 at the 
outer rational surface (ro = rs 2 ) saturates the DTM, as shown 
by contours of the helical flux ip* (left) and perpendicular 
flow (right) at the final simulation time of t = 750. 


differential rotation simulations.^ Once recoupled, the 
separatrix merging event proceeds with only minor devi¬ 
ation from the force-free system. Thus locating a drift at 
the inner, sub-dominant rational surface results in more 
system kinetic energy (due to plasma flows near the in¬ 
ner surface) but is not an effective means of slowing DTM 
mode growth. 


B. Outer rational surface 


Locating a strong diamagnetic drift at the outer reso¬ 
nant surface strongly stabilizes the nonlinear DTM. The 
early time growth of Em in Figure shows large oscil¬ 
lations, indicating that the two tearing surfaces are ini¬ 
tially decoupled. At later times these fluctuations con¬ 
tinue, but at a smaller amplitude compared to the total 
perturbed magnetic energy. 


Considering the state of the system at the last simu¬ 
lation time {t = 750, Figure 10), a significant amount 
of flux remains between the two q = 2 surfaces and the 
annular current ring is intact. The plasma flow is largely 
circulating inside the outer islands rather than between 
the two surfaces. The flattening of the energy growth, to¬ 
gether with the relaxed magnetic structure in Figure [TOl 
shows that the fundamental m = 2, n = 1 double-tearing 
mode is effectively saturated. Intermittent reconnection 
activity causes late time fluctuations of Ek and Em as the 
system oscillates, but does not significantly disturbed the 
saturated state. At these late times the source electric 
field causes slow growth of the perturbed perturbed mag¬ 
netic energy, as it drives the system back toward equilib¬ 
rium. By comparing to simulations without the source 
field we have determined that it has the effect of ‘pump¬ 
ing’ the simulation with energy rather than allowing it 
to fully relax. 




















t=250 t=390 




FIG. 11. Simulation times t = 250 (kinetic energy roll¬ 
over) and t = 390 (last simulation time) for an equal drift 
Ldt =0.1 DTM. Nonlinear enhancement of the pressure gra¬ 
dient severely inhibits reconnection so that flux remains be¬ 
tween the two rational surfaces, but cannot halt the structural 
instability. 


C. Equal drift 



The nonlinear development of the equal drift profile, 
as viewed through the energy evolution, is significantly 
more complicated than the previous case. growth 
slows significantly near t « 250 (Figu re [^ , then rises at 
a reduced, fluctuating rate. In Figure [Tl|( at this crest in 
Ek) a signihcant amount of flux remains between the two 
tearing surfaces. This roll-over does not, therefore, corre¬ 
spond to a separatrix merging event. Even at the end of 
our simulation {t = 390), when the inner current sheets 
approach each other across the magnetic axis, a small 
amount of flux remains between the two separatrices. 

The decrease in and E^ growth is a result of non¬ 
linear evolution of the pressure profile. Figure [T^ shows 
cuts of drp/p (the pressure and density contribution to 
Eqn. |14[ ) across the inner and outer current sheets at time 
t = 250 (the roll-over point) compared to the equilibrium 
profile. The nonlinear growth of the magnetic islands dis¬ 
torts the equilibrium pressure gradient and results in a 
significant enhancement to the diamagnetic drift within 
the tearing layers. As a consequence, reconnection is 
highly suppressed. 

Similar nonlinear enhancement of the pressure gradi¬ 
ent has been proposed a s a m echanism for saturation 
TO = 1 kink-tearing mode.l^Sin] jn our DTM simulation, 
however, this nonlinear enhancement of w* does not lead 
to complete stabilization. The magnetic islands con¬ 
tinue to evolve, and intermittent reconnection occurs. 
This continued growth of the kinetic and magnetic ener¬ 
gies (Fig. and deformation of the separatrices (Fig. [TT]) 
is a consequence of the large islands sizes required to en¬ 
hance the pressure gradient sufficiently to cut off recon¬ 
nection. The magnetic structure has become unstable, 
similar to pre vious simulations of explosive-type double¬ 
tearing modes,f2lIl3l instability continues to 

develop instead of saturating. In this respect the equal 
drift profile is less desirable than the outer drift, which 
stabilized the mode before significant deformation. 


FIG. 12. Cuts of drp/p across the inner and outer current 
sheets at t = 250 show significant enhancement of the pressure 
gradient compared to the equilibrium. 


VI. CONCLUSION 

In this work we have shown that electron diamagnetic 
drifts can have a stabilizing effect on the to = 2 , n = 1 
double-tearing mode. Their efficacy depends, however, 
on where they are located. Linearly we were only able to 
achieve a substantial decrease in the DTM growth rate 
by localizing a strong pressure gradient at the outer ra¬ 
tional surface. This class of prohle combined the decou¬ 
pling properties of differential plasma rotation with the 
local stabilizing benefits of the diamagnetic drift on the 
dominant tearing layer. The equal and inner drift pro¬ 
files contained only one of these stabilizing effects, and 
thus could not overcome the increased ideal MHD driving 
caused by the pressure gradient. 

Nonlinearly we found the decoupling and local stabi¬ 
lization effects of the outer drift profile resulted in satura¬ 
tion of the DTM at finite amplitude. The preservation of 
the annular current ring in this simulation indicates that 
diamagnetic drifts may act as a mechanism for stabilizing 
off-axis sawtooth crashes. We also found that nonlinear 
enhancement of the pressure gradient in an equal drift 
profile was able to significantly slow the growth of the 
instability. The large islands necessary to halt reconnec¬ 
tion, however, resulted in unstable magnetic structure. 

Our work shows that when hnite Larmor radius ef¬ 
fects are included the presence of DTM activity depends 
strongly on the details of the plasma pressure. The equi¬ 
libria considered in this work are highly constrained, and 
more realistic profiles will likely include some interme¬ 
diate form of our results. Whether such profiles are a 
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viable mechanism for avoiding DTM driven off-axis saw¬ 
tooth behavior in RMS devices is currently unclear. Ex¬ 
tending these results to high mode number DTMs and 
toroidal geometry would allow better comparison to ex¬ 
perimental data, as would consideration of non-constant 
temperatures and hot ions. We have provided strong ev¬ 
idence, however, for the dependence of DTM activity on 
the location of steep pressure profiles in the plasma. 
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